ELSEVIER 


Available  online  at  www.sciencedirect.com  - 

JOURNAL  OF 

ScienceDirect  POWER 

SOURCES 

Journal  of  Power  Sources  111  (2008)  500-505 

www.elsevier.com /locate /jpowsour 

Short  communication 


The  reaction  pathway  and  rate-limiting  step  of 
dehydrogenation  of  the  LiHN2  +  LiH  mixture 

Leon  L.  Shaw*,  William  Osborn, 

Tippawan  Markmaitree,  Xuefei  Wan 

Department  of  Chemical,  Materials  and  Biomolecular  Engineering, 

University  of  Connecticut,  Storrs,  CT,  USA 

Received  8  October  2007;  received  in  revised  form  8  November  2007;  accepted  9  November  2007 
Available  online  21  November  2007 


Abstract 

The  reaction  pathway  and  rate-limiting  step  of  dehydrogenation  of  the  LiNH2  +  LiH  mixture  have  been  investigated.  The  study  reveals  that 
dehydrogenation  of  the  LiNH2  +  LiH  mixture  is  diffusion-controlled  and  the  rate-limiting  step  is  NH3  diffusion  through  the  Li2NH  product  layer 
outside  the  LiNH2  shrinking  core.  This  phenomenon  is  explained  based  on  a  model  describing  the  major  steps  of  the  dehydriding  reaction  of 
the  mixture,  and  related  to  the  evidence  obtained  from  X-ray  diffraction  and  specific  surface  area  measurements  of  the  mixture  before  and  after 
isothermal  hydrogen  uptake/release  cycles  at  high  homologous  temperatures. 

©  2007  Elsevier  B.V.  All  rights  reserved. 

Keywords:  Hydrogen  storage  materials;  LiNH2;  LiH;  Reaction  kinetic;  Reaction  mechanism 


1.  Introduction 

Since  the  first  report  by  Chen  et  al.  [1]  in  2002,  studies  on 
hydrogen  sorption  and  desorption  behavior  and  mechanisms  of 
the  lithium  amide  (LiNH2)  and  lithium  hydride  (LiH)  mixture 
have  been  very  active  [2-15].  It  is  generally  agreed  that  the 
overall  dehydriding  reaction  of  this  system  is  expressed  as  [1]: 

LiNH2  +  LiH  =  Li2NH  +  H2  (1) 

This  reaction  can  theoretically  absorb  and  desorb  6.5wt.% 
hydrogen,  and  its  reaction  heat  has  been  calculated  to  be 
44.5  kJ mol-2  H2  [1],  although  a  recent  measurement  [9]  sug¬ 
gests  that  the  reaction  heat  might  be  65.6  kJ  mol-2  H2,  higher 
than  the  previous  theoretical  estimation  [1].  Additionally,  many 
studies  [4,5,11,16]  have  provided  strong  evidence  that  Reac¬ 
tion  (1)  proceeds  with  two  elementary  reactions.  First,  LiNH2 
decomposes,  as  shown  in  Reaction  (2): 

LiNH2  =  l/2Li2NH  +  1/2NH3  (2) 
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The  ammonia  (NH3)  from  Reaction  (2)  then  reacts  with  LiH  to 
form  LiNH2  again  and  liberate  H2,  as  shown  in  Reaction  (3): 

1/2NH3 -fl/2LiH  =  l/2LiNH2  +  l/2H2  (3) 

For  a  mixture  of  LiNH2  -1-LiH  (with  a  molar  ratio  of  1:1),  the 
reaction  would  continue  to  repeat  the  cycle  of  Reactions  (2)  and 
(3)  until  all  LiNH2  and  LiH  transform  to  lithium  imide  (Li2NH) 
and  H2  completely. 

It  has  been  demonstrated  that  Reaction  (3)  is  ultrafast  and 
can  take  place  in  the  order  of  microseconds  [4].  In  contrast, 
using  pure  LiNH2  as  the  starting  material,  it  has  been  shown 
that  Reaction  (2)  is  diffusion-controlled  and  proceeds  in  min¬ 
utes  [17].  It  is  interesting  to  note  that  both  Reactions  (2)  and 
(3)  have  solid  products,  yet  one  could  proceed  in  microseconds 
and  the  other  requires  much  longer  time  to  complete.  The  rea¬ 
son  for  such  a  large  difference  in  the  reaction  kinetics  between 
Reactions  (2)  and  (3)  has  not  been  studied.  The  reaction  kinet¬ 
ics  of  Reaction  (1),  i.e.,  the  combined  reaction  of  Reactions  (2) 
and  (3),  has  never  been  investigated  either.  This  study  has  been 
conducted  as  the  first  attempt  to  investigate  these  issues  and 
identify  the  rate-limiting  step  for  dehydrogenation  of  Reaction 
(1).  The  understanding  developed  from  this  study  will  shed  light 
on  how  to  increase  the  dehydrogenation  rate  of  the  LiNH2  +  LiH 
mixture. 
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2.  Materials  and  methods 

Lithium  amide  with  95%  purity  was  purchased  from  Alfa 
Aesar,  while  lithium  hydride  with  95%  purity  was  purchased 
from  Sigma-Aldrich.  The  LiNH2  and  LiH  mixture  was  prepared 
with  a  molar  ratio  of  1:1.1  according  to  Reaction  (1).  The  10% 
excess  of  LiH  was  added  to  minimize  the  loss  of  NH3  during  the 
dehydriding  process.  High-energy  ball  milling  was  conducted 
using  a  modified  Szegvari  attritor,  which  had  been  shown  pre¬ 
viously  to  be  effective  in  preventing  the  formation  of  the  dead 
zone  and  producing  uniform  milling  products  within  the  powder 
charge  [18].  Sample  loading  and  subsequent  ball  milling  were 
conducted  under  an  argon  atmosphere,  and  other  ball  milling 
conditions  can  be  found  elsewhere  [16]. 

After  ball  milling  for  3  h  at  room  temperature,  the 
LiNH2  +  LiH  mixture  was  then  subjected  to  various  character¬ 
izations  and  handled  in  a  glove-box  filled  with  Ar  of  99.999% 
purity.  The  characterization  included  (i)  the  specific  surface 
area  (SSA)  measurement  through  nitrogen  adsorption  at  77  K 
based  on  the  Brunauer-Emmett-Teller  (BET)  method  [19]  using 
a  gas  sorption  analyzer  (NOVA  1000),  (ii)  X-ray  diffraction 
(XRD)  using  a  D5005  ADAVANCE  diffractometer  with  Cu 
Ka  radiation,  and  (iii)  an  isothermal  soak/release  cycle  using 
a  Sieverts’-type  pressure-composition-isotherm  (PCI)  device 
(Advanced  Materials  Corporation).  Detailed  experimental  con¬ 
ditions  for  SSA  measurement  and  XRD  data  collection  can  be 
found  elsewhere  [16],  and  thus  will  not  be  repeated  here. 

The  isothermal  soak/release  cycle  entailed  a  1  h  sorption  at  a 
hydrogen  pressure  of  10  atm  and  a  subsequent  2.5  h  desorption 
under  an  evacuated  condition,  at  285  or  240  °C  and  repeated 
for  10  times.  The  sorption  segment  in  each  soak/release  cycle 
was  done  in  one  step,  i.e.,  exposing  the  sample  to  a  hydrogen 
pressure  of  10  atm  until  1  h  was  reached.  In  contrast,  the  des¬ 
orption  segment  in  each  soak/release  cycle  was  accomplished 
in  10  evacuation  steps  in  order  to  maintain  an  evacuated  condi¬ 
tion.  The  sample  chamber  was  closed  for  15  min  between  every 
two  consecutive  evacuation  steps  in  order  to  quantify  how  much 
hydrogen  was  released  in  the  next  evacuation  step.  During  each 
15  min  holding,  the  pressure  of  the  sample  chamber  increased 
gradually  because  of  desorption  of  hydrogen  from  the  sample 
and  was  recorded  to  calculate  the  hydrogen  desorbed  as  a  func¬ 
tion  of  time.  The  sum  of  the  hydrogen  desorbed  in  10  evacuation 
steps  represented  the  total  amount  of  hydrogen  desorbed  from 
the  sample  in  one  release  cycle  which  lasted  2.5  h  (i.e.,  15  min 
per  holding  x  10  steps).  The  pressure  of  the  sample  chamber  was 
less  than  0.003  atm  at  the  end  of  the  10th  holding  step.  Thus,  the 
total  amount  of  hydrogen  desorbed  in  one  release  cycle  could  be 
regarded  approximately  as  the  amount  of  hydrogen  released  in 
2.5  h  under  a  constant  temperature  (e.g.,  285  °C)  and  a  constant 
hydrogen  pressure  (i.e.,  0.003  atm). 

3.  Results  and  discussion 

Shown  in  Fig.  1  are  the  absorbed  and  desorbed  hydro¬ 
gen  (in  wt.%)  of  the  LiNH2  +  LiH  mixture  during  isothermal 
soak/release  cycles  at  285  °C.  It  is  clear  that  the  mixture  has 
a  rapid  hydrogen  sorption  rate  (i.e.,  approaching  the  theoreti- 


Time  (hr) 

Fig.  1.  Isothermal  hydrogen  uptake/release  cycles  of  the  LiNIL  +LiH  mixture 
at  285  °C:  (a)  the  overall  view  of  H2  absorbed  and  desorbed  as  a  function  of 
cycles,  and  (b)  the  close-up  of  the  5th  uptake/release  segment. 

cal  storage  capacity  in  ~5min)  and  slow  desorption  rate  (i.e., 
incomplete  release  of  hydrogen  in  2.5  h).  Because  of  its  slow  des¬ 
orption  rate,  the  mixture  can  only  desorb  ~4.4  wt.%  H2  in  each 
release  segment.  This  number  is  clearly  lower  than  its  theoretical 
storage  capacity  which  is  estimated  to  be  about  5.7  wt.%  H2  if  the 
extra  10  mol%  LiH  addition  and  the  presence  of  the  oxides  in  the 
starting  materials  are  considered.  As  a  result  of  the  un-released 
hydrogen  at  the  end  of  each  release  segment,  the  amount  of  the 
absorbed  hydrogen  can  only  be  4.4  wt.%  in  each  soak  segment. 
It  is  also  noted  that  little  degradation  in  the  kinetic  performance 
occurs  over  the  10  soak/release  cycles.  This  is  remarkable,  con¬ 
sidering  that  the  cyclic  temperature  (285  °C)  is  86%  of  LiN^’s 
melting  temperature  and  58%  of  LiH’ s  melting  temperature,  and 
that  the  holding  time  is  35  h.  It  is  well  known  that  holding  powder 
compacts  at  high  homologous  temperatures  (e.g.,  >0.5Tm  where 
Tm  is  the  melting  temperature  of  the  material)  for  several  hours 
is  sufficient  to  cause  annealing,  grain  growth,  and  even  densifi- 
cation  [20].  Thus,  the  unusual  cyclic  stability  in  the  soak/release 
performance  suggests  that  other  mechanism(s)  have  operated  to 
counter  balance  the  well-known  grain  growth  and  densification 
effects. 

Table  1  summarizes  the  crystallite  sizes  and  SSA  of  the 
LiNH2  +  LiH  mixture  before  and  after  isothermal  soak/release 
cycles.  Note  that  the  crystallite  sizes  of  LiH  and  LiNH2  after  ball 
milling,  determined  from  the  XRD  peak  broadening  through 
the  Scherrer  formula  [21]  with  the  correction  of  instrumental 
broadening,  are  in  nanoscales.  Thus,  it  can  be  concluded  that 
ball  milling  for  3  h  has  mixed  LiH  and  LiNH2  uniformly  at  the 
nanometer  scale.  It  is  further  noted  that  the  crystallite  size  of 


502 


L.L.  Shaw  et  al.  /  Journal  of  Power  Sources  177  (2008)  500-505 


Table  1 

XRD-determined  crystallite  sizes  of  LiH  and  LiNH2  and  the  measured  specific  surface  area  of  their  mixture 


Sample  conditions 

Crystallite  size 

Crystallite  size 

Specific  surface  area  of 

Equivalent  particle  size 

of  LiH  (nm) 

of  LiNH2  (nm) 

LiH  +  LiNH2  (m2g“2) 

of  LiH  +  LiNH2a  (nm) 

Before  isothermal  Cycling 

20.8 

13.3 

55 

51.8 

After  isothermal  Cycling 

20.0 

15.9 

10 

284.9 

a  The  equivalent  particle  size  of  the  (LiH  +  LiMb)  mixture  was  calculated  using  the  specific  surface  area  data  and  assuming  that  LiH  and  LiNH2  had  the  same 
particle  size. 


LiH  remains  almost  the  same  after  35  h  exposure  at  285  °C  and 
10  soak/release  cycles,  whereas  the  crystallite  size  of  LiNH2 
exhibits  only  a  slight  increase.  The  slight  or  no  growth  of  LiH 
and  LiNH2  nanocrystallites  is  truly  remarkable,  which  further 
confirms  that  mechanisms  opposing  grain  growth  have  been 
active  during  the  isothermal  soak/release  cycles  to  maintain 
the  nanocrystallites  of  LiH  and  LiNH2,  and  thus  offers  partial 
explanation  why  the  LiNH2  +  LiH  mixture  has  an  unusual  cyclic 
stability  in  the  soak/release  performance.  In  a  previous  study 
[22],  it  has  been  established  that  pure  LiH  particles  after  ball 
milling  with  the  identical  milling  condition  as  in  this  study  are 
agglomerates  of  single  crystals  with  the  median  crystallite  size  of 
12  nm,  which  grows  by  100%  after  the  ball-milled  LiH  particles 
expose  to  285  °C  for  1  h  in  an  argon  atmosphere.  Thus,  nano- 
LiH  particles  do  grow  as  expected,  when  exposed  to  285  °C 
without  soak/release  cycles  [22].  By  comparing  the  behav¬ 
iors  of  thermally  exposed  nano-LiH  particles  with  and  without 
soak/release  cycling,  it  can  be  concluded  that  soak/release 
cycles  do  have  effects  on  preventing  nano-LiH  crystallites  from 
growing. 

Note  that  Table  1  also  shows  a  large  decrease  in  the  specific 
surface  area  of  the  LiNH2  +  LiH  mixture  after  10  isothermal 
soak/release  cycles  at  285  °C.  The  decrease  in  the  SSA  can 
only  be  induced  by  (i)  crystal  growth,  (ii)  further  agglomeration 
of  multiple  particles,  or  (iii)  both.  However,  the  XRD  analy¬ 
sis  has  revealed  little  or  no  growth  of  nanocrystallites.  Thus, 
taking  the  SSA  and  crystallite  size  information  together,  it  can 
be  inferred  that  the  ball-milled  LiNH2  +  LiH  mixture  becomes 
further  agglomerated,  while  their  crystallite  sizes  remain  lit¬ 
tle  changed  after  exposing  to  isothermal  soak/release  cycles  at 
285  °C.  Additional  insights  can  be  attained  if  the  crystallite  sizes 
of  LiH  and  LiNH2  determined  from  XRD  are  compared  with 
the  equivalent  particle  size  of  the  LiNH2  +  LiH  mixture  calcu¬ 
lated  from  the  specific  surface  area  measurement.  As  shown  in 
Table  1 ,  the  equivalent  particle  size  of  the  LiNH2  +  LiH  mix¬ 
ture  is  larger  than  the  crystallite  sizes  of  both  LiH  and  LiNH2 
before  isothermal  cycling,  suggesting  the  presence  of  particle 
agglomeration.  This  result  is  in  excellent  agreement  with  pre¬ 
vious  transmission  electron  microscopy  (TEM)  analysis  [22], 
showing  agglomeration  of  single  crystals  of  pure  LiH  particles 
after  ball  milling  at  the  same  milling  condition  as  in  this  study. 
Table  1  also  reveals  that  agglomeration  becomes  worse  after 
isothermal  cycling  because  of  the  further  increased  equivalent 
particle  size  of  the  LiNH2  +  LiH  mixture  and  little  or  no  change 
in  the  crystallite  sizes  of  LiH  and  LiNH2.  However,  it  is  inter¬ 
esting  to  note  that  the  increased  agglomeration  does  not  result 
in  obvious  decreases  in  the  soak/release  performance  over  10 


(a)  LbNH 


Fig.  2.  Schematic  of  the  hydrogen  release  pathway  of  the  LiNH2  +  LiH  mixture: 
(a)  the  LbNH  product  from  Reaction  (2)  forms  a  continuous  shell  outside  the 
LiNH2  shrinking  core,  leading  to  a  reaction  rate  controlled  by  NH3  diffusion 
through  the  LbNH  product  layer,  and  (b)  the  LiMb  product  from  Reaction  (3) 
flakes  off  continuously,  resulting  in  direct  reaction  between  NH3  and  constantly 
regenerated  new  LiH  surface. 


cycles  in  35  h  (see  Fig.  1),  suggesting  that  the  surface  area  is  not 
the  rate-limiting  factor  in  the  soak/release  cycle  at  the  condition 
employed  in  this  study. 

In  order  to  identify  the  mechanism  responsible  for  maintain¬ 
ing  LiH  and  LiNH2  nanocrystallites  during  isothermal  cycling 
and  to  determine  the  rate-limiting  step  of  the  dehydriding  reac¬ 
tion  of  the  LiNH2  +  LiH  mixture,  the  desorption  rate  in  each 
release  segment  has  been  analyzed  using  various  kinetic  models 
that  are  possible  in  this  system.  As  established  in  Ref.  [17],  the 
rate  of  Reaction  (2)  is  controlled  by  diffusion  of  NH3  through  a 
porous  Li2NH  product  layer  (Fig.  2a).  Such  a  diffusion  is  neces¬ 
sary  for  continued  decomposition  of  LiNH2  at  the  LiNH2/Li2NH 
interface  and  driven  by  the  composition  gradient  of  high  NH3 
concentration  at  the  interface  and  low  NH3  concentration  at  the 
free  surface  of  L^NH  where  desorption  of  NH3  takes  place  [17]. 
Note  that  such  a  conclusion,  i.e.,  the  presence  of  a  porous  L^NH 
product  layer  on  a  shrinking  LiNH2  core  during  decomposition 
of  LiNH2,  is  consistent  with  the  expectation  if  the  changes  in 
the  volumes  of  the  solid  reactant  and  product  are  examined.  It 
is  well  known  that  the  Pilling-Bed worth  ratio  (P-B  ratio),  as 
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shown  in  Eq.  (4),  dictates  the  porosity  of  the  oxide  scale  and  its 
adhesion  to  the  metal  during  metal  oxidation  [23,24]. 

P  -  B  ratio  =  —  (4) 

V 

*  m 

where  V0  and  Vm  are  the  volumes  of  the  oxide  scale  and  the 
metal  oxidized,  respectively.  When  the  P-B  ratio  is  less  than 
1,  the  oxide  scale  tends  to  be  adherent  and  porous.  When  the 
P-B  ratio  is  between  1  and  2,  the  oxide  scale  is  adherent  and 
can  fully  cover  the  metal  surface  although  there  exist  compres¬ 
sive  stresses  in  the  scale.  When  the  P-B  ratio  is  large  than  2, 
the  oxide  scale  often  cracks  and  flakes  off  because  of  the  large 
compressive  stress  in  the  scale.  The  molar  volume  of  LiNH2  is 
19.44  cm3  mol-2,  whereas  the  corresponding  value  for  L^NH 
is  19.51  cm3  mol-2.  Based  on  the  formula  of  Reaction  (2)  and 
the  molar  volumes  of  LiNH2  and  L^NH,  the  P-B  ratio  for  Reac¬ 
tion  (2)  is  found  to  be  0.5,  thereby  allowing  the  formation  of  an 
adherent  yet  porous  L^NH  product  layer  through  which  NH3 
has  to  diffuse  and  then  desorb  from  the  L^NH  surface  in  order 
for  decomposition  of  LiNH2  to  continue.  When  the  P-B  ratio  of 
Reaction  (3)  is  considered  with  the  molar  volume  of  LiH  being 
9.70  cm3  mol-2  and  the  molar  volume  of  LiNH2  defined  above, 
it  is  found  that  the  P-B  ratio  is  2.0.  Thus,  the  LiNH2  product 
layer  from  Reaction  (3),  if  present  on  the  surface  of  the  LiH 
core,  would  have  extremely  large  compressive  stresses  and  tend 
to  flake  off  easily.  Therefore,  it  is  proposed  that  it  is  the  con¬ 
tinuous  flaking  off  of  the  LiNH2  product  layer  that  constantly 
provides  new  LiH  surface  for  reaction  with  NH3  and  thus  Reac¬ 
tion  (3)  can  proceed  in  the  order  of  microseconds,  as  reported 
in  Ref.  [4]. 

The  reaction  pathways  for  Reactions  (2)  and  (3)  proposed 
above  are  shown  schematically  in  Lig.  2.  Based  on  the  proposed 
reaction  pathways,  Reaction  (2)  can  at  least  be  described  with 
the  following  three  major  steps: 


LiNH2  -*  -Li2NH  +  -NH3# 

(5) 

1  1 

-nh3#  ->  -nh3* 

(6) 

1  1  1 

-nh3*  -*  -nh3  +  -* 

(7) 

Eq.  (5)  represents  decomposition  of  LiNH2  to  L^NH  and  NH3  at 
the  LiNH2/Li2NH  interface  with  the  NH3  absorbed  at  an  active 
interfacial  site,  #.  Eq.  (6)  stands  for  diffusion  of  NH3  through 
the  porous  L^NH  product  layer  from  an  active  interfacial  site, 
#,  to  an  active  L^NH  surface  site,  *  (see  Fig.  2).  Finally,  Eq. 
(7)  describes  desorption  of  NH3  from  the  L^NH  surface,  leav¬ 
ing  behind  an  active  surface  site.  Similarly,  Reaction  (3)  can  be 
described  with  at  least  three  major  steps: 


NH3  +  -  A  ->  -NH3A 

(8) 

1  1  1 

NH3A  +  -LiH  ->  -LiNH2  +  -H2A 

(9) 

Here  Eq.  (8)  represents  adsorption  of  NH3  on  an  active  LiH 
surface  site,  A,  whereas  Eq.  (9)  represents  the  absorbed  NH3 
reacting  with  LiH  to  form  LiNH2  flakes  and  H2  absorbed  at  an 
active  LiH  surface  site,  A.  Finally,  Eq.  (10)  describes  desorption 
of  H2  from  the  LiH  surface,  leaving  behind  an  active  surface  site. 

In  principle,  a  rate  equation  can  be  written  for  each  major 
step  listed  above.  However,  experiments  have  established  that 
Reaction  (3)  is  much  faster  than  Reaction  (2)  [4,17].  Thus,  it 
can  be  approximated  that  every  step  in  Reaction  (3)  is  in  equi¬ 
librium,  while  the  rate  of  Reaction  (1)  is  dictated  by  the  slow 
step(s)  of  Reaction  (2).  Therefore,  only  the  rate  of  each  step  in 
Reaction  (2)  needs  to  be  investigated  in  order  to  identify  the 
rate-limiting  step  of  Reaction  (1).  This  can  be  done  by  examin¬ 
ing  the  fraction  of  H2  released,/,  as  a  function  of  holding  time,  t , 
in  each  release  segment.  If  the  rate  of  Reaction  (1)  is  controlled 
by  Eq.  (5),  the  fraction  of  H2  released  as  a  function  of  hold¬ 
ing  time  can  be  derived  from  the  shrinking-core  model  (Fig.  2) 
with  a  LiNH2/Li2NH  interface  moving  at  a  constant  speed.  This 
relation  is  found  to  have  the  following  format  [25]: 

(l-/)1/3  =  l-^r  (11) 

However,  if  the  rate  of  Reaction  (1)  is  controlled  by  Eq.  (6), 
then  the  fraction  of  H2  released  will  follow  a  parabolic-rate  law 
derived  from  Fick’s  second  law  [26]  with  the  following  formula 
[17]: 

(1  -  /)1/3  =  1  -  ^-r1/2  (12) 

Finally,  if  the  rate  of  Reaction  (1)  is  controlled  by  Eq.  (7),  i.e., 
desorption  of  NH3  into  a  chamber  with  a  constant  pressure,  the 
fraction  of  H2  released  will  exhibit  a  linear  relationship  between 
/and  t  [27],  as  shown  by  Eq.  (13): 

/  =  kjRh  (13) 

The  £5,  &6  and  kj  above  are  rate  constants  for  Eqs.  (1 1),  (12)  and 
( 1 3),  respectively,  and  R  is  the  average  radius  of  LiNH2  particles. 

To  compare  the  rate  Eqs.  (11),  (12)  and  (13)  with  the  exper¬ 
imental  data  shown  in  Fig.  1,  these  rate  equations  need  to  be 
modified  to  take  into  account  the  incomplete  release  of  H2  at  the 
end  of  each  release  segment.  For  example,  Eq.  (11)  should  be 
modified  to 

(1  -  /7)1/3  =  1  -  (14) 

where  /  is  the  fraction  of  H2  released  at  the  end  of  the  release 
segment  in  reference  to  the  hydrogen  storage  capacity  of  the 
system,  which  is  5.7  wt.%  H2  as  discussed  previously,  whereas 
/is  the  fraction  of  H2  released  at  any  given  time  in  reference  to 
the  end  of  the  release  segment.  /=  1  at  the  end  of  each  release 
segment,  whereas  /  =  1  only  if  the  H2  released  at  the  end  of 
the  release  segment  equals  the  hydrogen  storage  capacity  of  the 
system.  Otherwise,  /  is  smaller  than  1.  Similarly,  Eqs.  (12)  and 
(13)  should  be  modified,  respectively,  to  be  Eqs.  (15)  and  (16). 
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Fig.  3.  The  amount  of  H2  released  in  5th  and  9th  release  segments  shown  in 
Fig.  1  is  plotted  in  (a)  the  (1  —  ff)m  vs.  tm  fashion,  as  defined  by  Eq.  (15),  (b) 
the  (1  —  fj)1/3  vs.  t  fashion,  as  defined  by  Eq.  (14),  and  (c)  the//vs.  t  fashion,  as 
defined  by  Eq.  (16).  Linear  equations  from  curve  fitting  and  the  corresponding 
R 2  values  are  indicated. 

f'f  =  hR2t  06) 

Fig.  3  compares  the  plots  of  the  experimental  data  of  the 
5th  and  9th  hydrogen  release  segments  presented  in  differ¬ 
ent  axes  defined  by  Eqs.  (14),  (15)  and  (16).  Clearly,  the 
experimental  data  fits  the  diffusion-controlled  model  (Fig.  3a) 
far  better  than  the  moving-boundary-controlled  (Fig.  3b)  and 
the  NH3 -desorption-controlled  model  (Fig.  3c)  because  the 
diffusion-controlled  model  has  the  highest  R 2  value  obtained 
from  curve  fitting  through  the  least-squares  method.  This  con¬ 
clusion  applies  to  every  release  segment  in  Fig.  1  with  the  5th 


Fig.  4.  A  comparison  of  the  amounts  of  H2  released  as  a  function  of  time  in  the 
5th  release  segments  of  the  isothermal  soak/release  cycles  at  240  and  285  °C, 
plotted  in  the  (1  —  ff)1/3  vs.  tm  fashion,  as  defined  by  Eq.  (15). 


and  9th  segments  shown  in  Fig.  3.  In  spite  of  its  highest  R 2  value 
(Fig.  3a),  the  deviation  of  the  experimental  data  from  the  linear 
relationship  between  (1  —ff)113  and  tm  is  present.  This  devia¬ 
tion  is  due  to  the  gradual  increase  in  the  pressure  of  the  sample 
chamber  during  the  15  min  holding  period  in  each  evacuation 
step,  especially  at  the  early  stage  of  the  release  segment  when 
the  diffusion  distance  is  very  short  and  the  pressure  rise  is  rel¬ 
atively  quick.  However,  when  the  isothermal  soak/release  cycle 
is  performed  at  lower  temperatures  and  thus  both  the  pressure 
rise  and  the  diffusion  rate  become  slower,  the  linearity  improves 
substantially.  Fig.  4  shows  such  a  situation,  revealing  that  the 
R2  value  almost  becomes  1 .0  when  the  isothermal  soak/release 
cycle  is  conducted  at  240  °C.  Thus,  this  analysis  unambiguously 
reveals  that  the  rate-limiting  step  of  Reaction  (1)  is  diffusion  of 
NH3  through  the  porous  L^NH  product  layer. 

4.  Concluding  remarks 

The  present  set  of  experiments  reveals  that  nanocrystalline 
LiH  and  ONH2  particles  generated  via  high-energy  ball  milling 
exhibit  little  or  no  growth  in  nanocrystallite  sizes  after  10  hydro¬ 
gen  soak/release  cycles  at  285  °C  for  35  h.  However,  the  surface 
area  of  the  ONH2  +  LiH  powder  mixture  does  decrease,  indi¬ 
cating  an  increase  in  particle  agglomeration  during  isothermal 
cycling.  Nevertheless,  the  isothermal  hydrogen  soak/release  per¬ 
formance  exhibits  little  or  no  degradation.  The  unusual  cyclic 
stability  of  the  nanocrystalline  LiNH2  +  LiH  mixture  has  been 
attributed  to  (i)  repeated  nucleation  of  solid  phases  during 
hydriding  and  dehydriding  reactions  and  (ii)  the  continuous 
flaking  off  of  the  LiNH2  product  layer  from  the  surface  of  the 
LiH  reactant.  The  latter,  in  conjunction  with  the  presence  of  H2 
gas,  effectively  prevents  substantial  decrease  in  the  specific  sur¬ 
face  area  and  densification  during  isothermal  cycling  for  35  h 
at  very  high  homologous  temperatures  (i.e.,  86%  of  LfiNTb’s 
melting  temperature  and  58%  of  LiH’s  melting  temperature). 
Dehydrogenation  of  the  LiNH2  +  LiH  mixture  is  found  to  be 
diffusion-controlled  and  the  rate-limiting  step  is  identified  to  be 
diffusion  of  NH3  through  the  porous  L^NH  product  layer  out- 
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side  the  ONH2  shrinking  core.  A  reaction  pathway  has  been 
proposed  to  explain  ultrafast  reaction  between  LiH  and  NH3  as 
well  as  slow  decomposition  of  LiNH2  to  L^NH  and  NH3 .  The 
model  proposed  is  in  good  agreement  with  experimental  obser¬ 
vations  including  (i)  the  change  in  the  volume  of  the  solid  before 
and  after  reactions,  (ii)  the  limited  increase  in  particle  agglom¬ 
eration  over  35  h  at  very  high  homologous  temperatures,  and 
(iii)  the  substantially  different  reaction  rates  between  Reactions 
(2)  and  (3)  although  both  reactions  generate  solid  and  gaseous 
products. 
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